Coarse grained force field for the molecular simulation of natural gases and condensates

AbstractThe atomistically-detailed molecular modelling of petroleum fluids is challenging, amongst other aspects, due to the very diverse multicomponent and asymmetric nature of the mixtures in question. Complicating matters further, the time scales for many important processes can be much larger than the current and foreseeable capacity of modern computers running fully-atomistic models. To overcome these limitations, a coarse grained (CG) model is proposed where some of the less-important degrees of freedom are safely integrated out, leaving as key parameters the average energy levels, the molecular conformations and the range of the Mie intermolecular potentials employed as the basis of the model. The parametrization is performed by using an analytical equation of state of the statistical associating fluid theory (SAFT) family to link the potential parameters to macroscopically observed thermophysical properties. The parameters found through this top-down approach are used directly in molecular dynamics simulations of multi-component multi-phase systems. The procedure is exemplified by calculating the phase envelope of the methane–decane binary and of two synthetic light condensate mixtures. A methodology based on the discrete expansion of a mixture is used to determine the bubble points of these latter mixtures, with an excellent agreement to experimental data. The model presented is entirely predictive and an abridged table of parameters for some fluids of interest is provided

Formation of Methane Hydrate in the Presence of Natural and Synthetic Nanoparticles

Natural gas hydrates occur widely on the ocean-bed and in permafrost regions, and have potential as an untapped energy resource. Their formation and growth, however, poses major problems for the energy sector due to their tendency to block oil and gas pipelines, whereas their melting is viewed as a potential contributor to climate change. Although recent advances have been made in understanding bulk methane hydrate formation, the effect of impurity particles, which are always present under conditions relevant to industry and the environment, remains an open question. Here we present results from neutron scattering experiments and molecular dynamics simulations that show that the formation of methane hydrate is insensitive to the addition of a wide range of impurity particles. Our analysis shows that this is due to the different chemical natures of methane and water, with methane generally excluded from the volume surrounding the nanoparticles. This has important consequences for our understanding of the mechanism of hydrate nucleation and the design of new inhibitor molecules

Meso- and macrozooplankton communities in the Weddell Sea, Antarctica

The present paper describes composition and abundance of meso- and macrozooplankton in the epipelagic zone of the Weddell Sea and gives a systematic review of encountered species regarding results of earlier expeditions. Material was sampled from 6 February to 10 March 1983 from RV Polarstern with a RMT 1+8 m (320 and 4500 μm mesh size). In agreement with topography and water mass distribution three distinct communities were defined, clearly separated by cluster analysis: The Southern Shelf Community has lowest abundances (approx. 9000 ind./1000 m3). Euphausia crystallorophias and Metridia gerlachei are predominating. Compared with the low overall abundance the number of regularly occurring species is high (55) due to many neritic forms. Herbivores and omnivores are dominating (58% and 35%). The North-eastern Shelf Community has highest abundances (about 31 000 ind./1000 m3). It is predominated by copepodites I–III of Calanus propinquus and Calanoides acutus (61%). The faunal composition is characterized by both oceanic and neritic species (64). Fine-filter feeders are prevailing (65%). The Oceanic Community has a mean abundance of approximately 23 000 ind./1000 m3, consisting of 61 species. Dominances are not as pronounced as in the shelf communities. Apart from abundant species like Calanus propinquus, Calanoides acutus, Metridia gerlachei, Oithona spp. and Oncaea spp. many typical inhabitants of the Eastwind Drift are encountered. All feeding types have about the same importance in the Oceanic Community

Aggregation behavior of model asphaltenes revealed from large-scale coarse-grained molecular simulations

Fully atomistic simulations of models of asphaltenes in simple solvents have allowed the study of trends in aggregation phenomena and the understanding of the role that molecular structure plays therein. However, the detail included at this scale of molecular modeling is at odds with the required spatial and temporal resolution needed to fully understand the asphaltene aggregation. The computational cost required to explore the relevant scales can be reduced by employing coarse-grained (CG) models, which consist of lumping a few atoms into a single segment that is characterised by effective interac- tions. In this work CG force fields developed via the SAFT-γ [Müller, E.A., Jackson, G. (2014) Annu. Rev. Chem. Biomolec. Eng., 5, 405–427] equation of state (EoS) provide a reliable pathway to link the molecular description with macroscopic thermophysical data. A recent modification of the SAFT-VR EoS [Müller, E.A. and Mejía, A. (2017) Langmuir, 33, 11518–11529], that allows parametrizing homonuclear rings, is selected as the starting point to propose CG models for polycyclic aromatic hydrocarbons (PAHs). The new aromatic-core parameters, along with others published for simpler organic molecules, are adopted for the construction of asphaltene models by combining different chemical moieties in a group-contribution fashion. We apply the procedure to two previously reported asphaltene models and perform Molecular Dynamics simulations to validate the coarse-grained representation against benchmark systems of 27 asphaltenes in pure solvent (toluene or heptane) described in a fully atomistic fashion. An excellent match between both levels of description is observed for cluster size, radii of gyration, and relative-shape-anisotropy-factor distributions. We exploit the advantages of the CG representation by simulating systems containing up to 2000 asphaltene molecules in explicit solvent investigating the effect of asphaltene concentration, solvent composi- tion, and temperature on aggregation. Upon employing large systems facilitated by the CG models, we observe stable continuous distributions of molecular aggregates at con- ditions away from the two-phase precipitation point. As a further example application, a widely accepted interpretation of cluster-size distributions in asphaltenic systems is challenged by performing system-size tests, reversibility proofs and time-dependence analysis. The coarse-graining procedure proposed is seen to be general and predictive, hence, can be applied to other asphaltenic molecular structures

Multiscale Approach Linking Self-Aggregation and Surface Interactions of Synthesized Foulants to Fouling Mitigation Strategies

© 2019 American Chemical Society. Fouling of oil-exposed surfaces remains a crucial issue as a result of the continued importance of oil as the world's primary energy source. The key perpetrators in crude oil fouling have been identified as asphaltenes, a poorly described mixture of diverse polyfunctional molecules that form part of the heaviest fractions of oil. Asphaltenes are responsible for a decrease in oil production and energy efficiency and an increase in the risk of environmental hazards. Hence, understanding and managing systems that are prone to fouling is of great value but constitutes a challenge as a result of their complexity. In an effort to reduce that complexity, a study of a synthesized foulant of archipelago structure is presented. A critical perspective on previously described solubility and aggregation mechanisms (e.g., critical nanoaggrerate concentration and critical clustering concentration) is offered because the characterized system favors a continuous distribution of n-mers instead. A battery of experimental and modeling techniques have been employed to link the bulk and interfacial behavior of a representative foulant monomer to effective fouling mitigation strategies. This systematic approach defines a new multiscale methodology in the investigation of fouling systems